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Abstract 

Already in an early disease stage, patients with chronic obstructive pulmonary disease (COPD) are confronted with impaired 
skeletal muscle function and physical performance due to a loss of oxidative type I muscle fibers and oxidative capacity (i.e. 
oxidative phenotype; Oxphen). Physical activity is a well-known stimulus of muscle Oxphen and crucial for its maintenance. 
We hypothesized that a blunted response of Oxphen genes to an acute bout of exercise could contribute to decreased 
Oxphen in COPD. For this, 28 patients with less advanced COPD (age 65±7 yrs, FEV, 59±16% predicted) and 15 age- and 
gender-matched healthy controls performed an incremental cycle ergometry test. The Oxphen response to exercise was 
determined by the measurement of gene expression levels of Oxphen markers in pre and 4h-post exercise quadriceps 
biopsies. Because exercise-induced hypoxia and oxidative stress may interfere with Oxphen response, oxygen saturation 
and oxidative stress markers were assessed as well. Regardless of oxygen desaturation and absolute exercise intensities, the 
Oxphen regulatory response to exercise was comparable between COPD patients and controls with no evidence of 
increased oxidative stress. In conclusion, the muscle Oxphen regulatory response to acute exercise is not blunted in less 
advanced COPD, regardless of exercise-induced hypoxia. Hence, this study provides further rationale for incorporation of 
exercise training as integrated part of disease management to prevent or slow down loss of muscle Oxphen and related 
functional impairment in COPD. 
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Introduction 

Loss of skeletal muscle oxidative phenotype (Oxphen) is 
prevalent in chronic obstructive pulmonary disease (COPD) [1- 
3]. It includes a proportional shift from slow-oxidative type I 
muscle fibers towards the fast-glycolytic type II fibers, associated 
with a reduced capacity of oxidative metabolism and in advanced 
disease also mitochondrial dysfunction [4]. Loss of muscle Oxphen 
is related to functional impairments, such as a reduction in muscle 
endurance [4], whole body exercise capacity, and mechanical 
efficiency [5], and it has been proposed as a driver of 
cardiovascular and metabolic risk [6] and cachexia [7]. Loss of 
muscle Oxphen is most prominent in severe COPD, but we 
recently showed that the process is already ongoing in patients 
with less advanced COPD [8]. Hence, this earlier stage of the 



disease is of particular importance to gain mechanistic insight into 
the onset of Oxphen loss in COPD. 

Exercise is an important stimulus for maintenance of muscle 
Oxphen. Because patients with COPD are physically less active 
[9], it is tempting to attribute Oxphen loss in these patients to 
disuse. However, several recent reports could not show such an 
association between physical activity level and Oxphen [8,10,1 1]. 
Physiologically, muscle Oxphen is progressively stimulated and 
maintained by repeated bouts of exercise through induction of 
among others peroxisome proliferator-activated receptor (PPAR) 
coactivator- 1 a (PGC-lot), PPAR-a and -8, and mitochondrial 
transcription factor A (TFAM) [12,13]. We hypothesized that the 
response of Oxphen regulation to acute bouts of exercise is blunted 
in COPD, which could contribute to a loss of muscle Oxphen 
irrespective of physical activity level. To test this hypothesis we 
compared Oxphen markers and their regulators in muscle biopsies 
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obtained before and after an acute bout of exercise, between less 
advanced COPD patients and healthy controls who were matched 
for age and body mass index (BMI). 

Hypoxia is known to increase muscle glycolytic capacity relative 
to oxidative capacity [14,15]. In humans, hypoxic exposure also 
induces loss of muscle oxidative capacity by reducing mitochon- 
drial volume (reviewed in [16]). Exercise-induced oxygen desat- 
uration even occurs in some patients with less advanced COPD, 
potentially reflecting or leading to hypoxia in the exercising 
muscle. Experimental research in mice has shown that exposure to 
acute hypoxia can also induce oxidative stress in skeletal muscle 
mitochondria [17]. Indeed, elevated levels of oxidative stress 
markers have been shown in muscles of COPD patients with 
chronic hypoxemia [18]. Moreover, enhanced exercise-induced 
oxidative stress has frequently been observed in muscle of patients 
with severe COPD [19], particularly in the wasted subgroup [20]. 
Interestingly, in severe patients, oxidative stress is inversely 
associated with muscle endurance [21] and in vitro permeabilized 
muscle fibers of patients with severe COPD produce significandy 
more reactive oxygen species than those of healthy controls [22]. 
This is possibly related to a fiber-type shift, since studies in rats 
have shown that mitochondria within type II fibers release 
significandy more H 2 0 2 than those within oxidative fibers [23]. 
As such, oxygen desaturation and/or muscle oxidative stress could 
mediate an anticipated blunted Oxphen response in COPD. A 
secondary objective was therefore to determine the involvement of 
exercise-induced oxygen desaturation and oxidative stress in the 
response to an acute exercise bout. 

Methods 

Ethics statement 

Written informed consent was obtained from all subjects and 
the ethical review board of the Maastricht University Medical 
Center+ approved the study (08-2-059). The study was registered 
at www.trialregister.nl as NTR1402. 

Patient characteristics 

The study population included 28 clinically stable mild to 
moderate COPD patients (FEV, 59±16% predicted; P a ,0 2 
9.3±1.0kPa) and 15 healthy controls. Recruitment of patients 
and controls and exclusion criteria were previously described [8]. 
Briefly, patients were recruited from the outpatient clinic of 
Maastricht University Medical Center+ (MUMC+, Maastricht, 
The Netherlands) and via advertisements in local newspapers. 
Exclusion criteria were long-term oxygen therapy, oral corticoste- 
roid use, acute exacerbation requiring hospital admission in the 
past eight weeks, rehabilitation in the past six months and known 
co-morbidities potentially interfering with study outcome param- 
eters. Healthy controls were recruited via advertising in local 
newspapers and absence of co-morbidities and airflow limitation 
was verified through history-taking by a physician, and pulmonary 
function tests. Healthy controls were matched to COPD patients 
on age, BMI and sex distribution. Study subjects were carefully 
characterized regarding body composition, lung function, smoking 
status and daily physical activity level, as reported previously [8]. 
The protocol used to determine quadriceps function on a 
dynamometer was previously described [8]. Briefly, under strong 
encouragement, subjects performed thirty sequential volitional 
maximal contractions at an angular velocity of 90°/ s. Isokinetic 
quadriceps endurance was determined as proportional decline in 
peak torques (relative to the highest peak torque) per repetition 
(%/rep). 



Experimental protocol 

To study the response of Oxphen regulators and oxidative stress 
to acute exercise, subjects performed an incremental load cycle 
ergometry test. Quadriceps muscle biopsies were obtained directly 
before and 4 h after the exercise test for Oxphen and oxidative 
stress analyses to obtain a maximal response of metabolic and 
oxidative stress-related genes after exercise. During the entire 
exercise test, oxygen saturation was measured by means of 
traditional finger pulse oximetry, as well as using a portable 
continuous-wave near-infrared spectroscopy (NIRS) system to 
examine oxygen saturation locally at the muscular level. Arterial 
punctures were performed before exercise and at maximal 
exertion for gas and lactate analyses. 

Acute exercise protocol 

All subjects performed an incremental load cycle ergometry test. 
The rationale for this protocol was to invoke a rapid exercise 
response in terms of Oxphen regulation particularly in the 
quadriceps muscle (see below). Peak oxygen uptake (I / 'o 2pt . ak ) 
and peak work load were determined as previously described [24] . 
Anaerobic threshold was determined from ergospirometry data 
using the v-slope method [25]. Arterial punctures of the radial 
artery at rest and at Fo 2pi , ak were available from n= 19 COPD 
patients and n = 1 1 healthy subjects. Arterial blood gas and lactate 
concentrations were analyzed on a Chiron blood gas analyzer 865 
(Chiron Diagnostics, Emeryville, CA). 

Arterial and peripheral muscle oxygenation 

Arterial oxygen saturation was continuously monitored using a 
finger pulse oximeter. By post hoc stratification the patients were 
divided into those who desaturated during exercise (•S p ,o 2 fall & 
4%; COPD D , n=13) and those who did not (5p,o 2 fall<4%; 
COPD ND , n = 15) according to ATS criteria [26]. Quadriceps 
muscle oxygenation was measured using the PortaMon portable 
continuous-wave near-infrared spectrophotometer (NIRS, Artinis 
Medical Systems, Zetten, The Netherlands). The probe was 
attached with Velcro straps on the vastus lateralis of the non- 
biopsied leg approximately 10 cm proximally from the knee joint, 
and covered with a black cloth to eliminate environmental light. 
Measurements started 3 minutes before exercise while the subject 
was still at rest and continued for 2-5 minutes after the exercise 
test until the signal had stabilized. Baseline values were established 
during 3 minutes of unloaded cycling before the incremental 
exercise test. Tissue oxygenation was measured continuously. 
Traces were analyzed with OxySoft Software (v2.1.2; Artinis 
Medical Systems). The tissue saturation index (TSI) was calculated 
using spatial resolved spectroscopy using the equation below, with 
[0 2 Hb] and [HHb] representing concentrations of oxygenated 
and deoxygenated hemoglobin/myoglobin. 



TSI = 100 x 



[O z Hb] 



[0 2 Hb] + [HHb] 



Values were reported as change from baseline (A). Using mono- 
exponential curve fitting, the time constant of muscle reoxygena- 
tion (tau, x) was determined from the slope of the TSI curve during 
the reoxygenation phase. Twelve subjects (4 COPD ND , 4 COPD D 
patients, and 4 healthy controls) were excluded from further NIRS 
analyses mainly because their skin fold at the site of the NIRS 
measurement was greater than 3 cm (which is the operational limit 
for the NIRS probe), or because the signal was poor or did not 
show a clear drop in TSI. 
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Muscle biopsy and analyses 

Biopsies were obtained at rest and 4 h after the exercise test 
from quadriceps muscle [vastus lateralis) of the dominant leg from 
two separate incisions approximately 2 cm apart using the needle 
biopsy technique [27]. The interval between the exercise test and 
the second biopsy was chosen after careful study of the literature to 
obtain a maximal response of metabolic genes after exercise (e.g., 
[28,29]). Muscle tissue was frozen in melting isopentane precooled 
in liquid nitrogen and stored in aliquots at — 80°C for histology 
and molecular analyses. 

Techniques for sample preparation for analyses of fiber-type 
composition, gene expression, oxidative stress markers, metabolic 
enzyme activity as well as the corresponding analyses have been 
previously described [8]. Briefly, fiber-type composition was 
determined on 5 u.m serial cryosections using antibodies against 
myosin heavy chain (MyHC) I, MyHC Ha (both from Develop- 
mental Studies Hybridoma Bank, University of Iowa, USA), and 
laminin (Sigma, Zwijndrecht, the Netherlands). Fiber-type classi- 
fication was aided by myosin ATPase-activity staining with acidic 
pre-incubation at pH 4.40 [30]. For qRT-PCR, RNA was 
extracted from 10-30 mg muscle tissue using the ToTALLY 
RNA™ Kit (Ambion Ltd., Foster City, CA, USA) according to 
the supplier's protocol, followed by genomic DNA removal with 
the RNeasy Mini Kit with RNase-free DNase (Qiagen, Venlo, 
The Netherlands). 400 ng RNA was reverse transcribed to cDNA 
with anchored oligo(dT) primers according to the supplier's 
protocol (Transcriptor First Strand cDNA Synthesis kit, Roche 
Diagnostics, Woerden, The Netherlands). qRT-PCR primers were 
designed based on Ensembl transcript sequences or selected from 
literature [31,32] and ordered from Sigma Genosys (Zwijndrecht, 
the Netherlands), with primer details shown in Table SI and Table 
S2 in File SI. Using these analyses, regulators and markers of 
oxidative and glycolytic metabolism as well as oxidative stress were 
assessed. Contents of antioxidant enzymes were determined with 
the western blot protocol using anti-catalase antibody (Calbio- 
chem, San Diego, CA, USA) and anti-manganese superoxide 
dismutase and anti-Cu/Zn superoxide dismutase antibodies (Santa 
Cruz Biotechnology, CA, USA). Protein content of oxidative stress 
markers were identified using antibodies for protein carbonylation 
(anti-2,4-DNP moiety antibody, Oxyblot kit, Chemicon Interna- 
tional Inc., Temecula, CA, USA), total protein nitration (anti-3- 
nitrotyrosine antibody, Invitrogen, Eugene, OR, USA), and total 
malondialdehyde (MDA)-protein adducts (anti-MDA antibody, 
Academy Biomedical Company Inc., Houston, TX, USA). For 
protein carbonylation, carbonyl groups in the protein side chains 
were first derivatized to 2,4-dinitrophenylhydrazone (DNP) using 
the Oxyblot kit (Chemicon International Inc., Temecula, CA, 
USA) according to the manufacturer's instructions. Exercise 
responses of gene expression and oxidative stress product variables 
were summarized with fold inductions, which were calculated by 
dividing post-exercise values by baseline values. Contents of 
representative subunits of oxidative phosphorylation (OXPHOS) 
complexes were determined by western blot using an anti-total 
OXPHOS antibody cocktail (MS604-300, Abeam, Cambridge, 
UK) with GAPDH (anti-GAPDH, #2118, Cell Signaling Tech- 
nology, Leiden, The Netherlands) as a loading control. Metabolic 
activities for citrate synthase (CS, EC 2.3.3.1) and 3-hydroxyacyl- 
CoA dehydrogenase (HADH, EC 1.1.1.35) were assayed spectro- 
photometrically (Multiskan Spectrum; Thermo Labsystems, 
Breda, The Netherlands) as previously described [33-35]. Absolute 
CS and HADH activities were normalized to total protein. 



Statistics 

The assumptions of normality and homogeneity of variances 
were checked for all experimental groups with the Shapiro-Wilk 
test and Levene's test. Three different sets of statistical compar- 
isons were made. First, baseline comparisons between the COPD 
patients and controls were tested using independent Student's t test 
or Mann-Whitney test as appropriate. Comparisons between 
COPD ND , COPD D and control groups were made using one-way 
ANOVA (post hoc Gabriel), Welch's ANOVA (post hoc Games- 
Howell) or Kruskal-Wallis non-parametric ANOVA (post hoc 
Mann-Whitney U with Bonferroni correction) as appropriate. 
Discontinuous variables were tested using Fisher's Exact test. 
These same tests were subsequently employed to analyze 
differences in exercise-induced response of skeletal muscle 
metabolic genes and oxidative stress markers between the groups 
(expressed as fold inductions), whereas paired Student's / test or 
Wilcoxon signed-rank test were performed to analyze such 
changes within each study group. Correlations were tested using 
Pearson's correlation coefficient, or Spearman's p in case of non- 
normally distributed data. Analyses were performed using IBM 
SPSS Statistics 20.0 (IBM Corp., Armonk, NY). A />-value<().()5 
was considered statistically significant. 

Results 

Pulmonary and anthropometric characteristics 

Characteristics of the study subjects are summarized in Table 1. 
COPD patients had moderate airflow obstruction (FEVi 59 ±16). 
Compared to healthy controls, the COPD patients had a lower 
resting arterial oxygen saturation (S p ,02 96 (95, 97) vs 98 (98, 
100)%) and daily physical activity level (214 (129, 300) vs 328 (276, 
489) counts/ min). COPD D patients did not differ from COPD ND 
patients in terms of pulmonary function indices, but had a lower 
BMI. 

Oxidative phenotype in quadriceps muscle 

At baseline, as reported previously [8] and listed in Table 2, 
Oxphen markers were lower in COPD patients than in controls, 
which was reflected by lower type I and higher type II fiber 
proportions, and reduced content of some subunits of the oxidative 
phosphorylation complexes. Oxphen markers were not different 
between patients who did or did not desaturate during exercise 
(Table 2). Moreover, quadriceps endurance was reduced in this 
COPD cohort (decline in peak torque -1.74 (-2.04, -1.28) in 
COPD vs - 1.09 (-1.57, -0.91)%/repetition in healthy controls), 
but it was not different between COPD D and COPD ND patients 
(-1.86 (-2.02, -1.35) in COPD D vs -1.61 (-2.14, -1.25)%/ 
repetition in COPD ND ). 

Cycle ergometry performance 

Data from the cycle ergometry test are listed in Table 3. As 
expected, peak work rate was lower in COPD patients than in 
controls (72 (56, 96) vs 156 (134, 255) Watt). No differences were 
found between COPD D and COPD ND patients. 

Exercise-induced metabolic gene expression 

The acute exercise clearly evoked an upregulation of genes 
known to be exercise-responsive, including PGC- la-related co- 
activator PRC, glycolysis-involved hexokinase II and the master 
regulator of the hypoxic response, HIF-lot (Figure 1). The COPD 
patients did not respond differently to exercise than controls 
regarding expression of genes associated with muscle oxidative 
metabolism, nor of genes involved in glycolytic metabolism or 
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Table 1. Main characteristics of subjects. 





Healthy controls 


COPD 


COPD ND 


COPD D 


N (M/F) 


15 (9/6) 


28 (16/12) 


15 (9/6) 


13 (7/6) 


Age, y 


65±6 


65±7 


66±7 


65±6 


Height, cm 


170±11 


169±9 


168±7 


169±11 


Weight, kg 


72±12 


71 ± 1 0 


74±10 


67±10 


BMI, kg/m 


24.9±3.3 


25.1 ±2.8 


26.4±2.5 


23.5±2.3 # 


FFMI, kg/m 


18.0±1.9 


17.5±1.8 


18.0±2.0 


16.9±1.5 


Smoking status current/former/never 


1/7/7 


10/18/0™ 


7/8/0" 


3/1 0/0* 


FEV,, % predicted 


113±15 


59±16™ 


61 ±13*™ 


56±19*" 


FVC, % predicted 


120±17 


104±22" 


98 ±20* 


110±24 


D^co- % predicted 


95±19 


51 ±16*" 


55 ±16*™ 


47±16*" 


5p,o 2 at rest, % 


98 (98, 100) 


96 (95, 97)"" 


97 (96, 97)*** 


95 (94, 98)*™ 


P a ,o 2 at rest, kPa 


12.1 ±1.1 


9.3 ±1.0*" 


9.5 ±0.8*** 


9.1 ±1.2*** 


P a ,co 2 at rest, kPa 


5.3 ±0.2 


5.2 ±0.5 


5.1 ±0.5 


5.2 ±0.4 


Daily physical activity, counts/min 


328 (276, 489) 


214 (129, 300)** 


208 (140, 278)" 


220 (108, 366) 



Values are expressed as mean±SD or median (25th percentile, 75th percentile). Abbreviations: BMI, body mass index; FFMI, fat-free mass index; FEVt, forced expiratory 
volume in one second; FVC, forced vital capacity; D LiC o, diffusion capacity of the lungs for carbon monoxide; Sp,o 2 , oxygen saturation measured via pulse oximetry. 
Significance of difference compared to controls: *p<0.05, **p<0.01, ™*p^0.001. Level of significance COPD D vs. COPD ND patients: # p<0.05. 
doi:1 0.1 371 /journal.pone.00901 50.t001 



hypoxia signaling. We also did not observe a different response in 
COPD D patients compared to COPD ND patients. 

Exercise response of oxidative stress markers 

At baseline, oxidative stress-induced protein modification levels 
such as carbonylation, tyrosine nitration and MDA adducts were 
not different between COPD patients and controls [8], nor 
between COPD D and COPD ND patients (data not shown). After 



exercise, gene expression levels of heme oxygenase-1 (HO-1) and 
manganese superoxide dismutase (MnSOD) were increased in 
both COPD patients and healthy controls (Figure 2A). An exercise 
response of the antioxidant enzymes was not detected at the 
protein level as manganese and copper/zinc superoxide dismutase 
and catalase protein content remained unchanged (Figure 2B). No 
differential response was shown between COPD D and COPD ND 
patients. 



Table 2. Main characteristics of vastus lateralis muscle fibers. 





Healthy controls 


COPD 


COPD ND 


COPD D 


Type 1 proportion, % 


62.7 (55.3, 68.3) 


39.7 (35.3, 48.6)*" 


40.3 (35.5, 62.3)* 


39.7 (29.8, 47.4)*** 


Type l/lla proportion, % 


8.1 (4.6, 11.0) 


5.6 (2.2, 9.5) 


5.3 (3.5, 9.0) 


6.4 (1.2, 12.7) 


Type ila proportion, % 


25.5 (16.8, 32.4) 


40.5 (25.9, 48.9)" 


32.0 (24.5, 46.1)* 


45.4 (39.7, 55.2)"** 


Type lla/llx proportion, % 


1.0 (0.3, 3.4) 


3.6 (1.1, 8.1) 


4.3 (1.0, 9.4) 


3.6 (0.7, 7.4) 


Type llx proportion, % 


0.5 (0.0, 2.7) 


2.3 (0.5, 9.4)* 


2.2 (0.2, 8.8) 


3.2 (0.6, 1 1 .0) 


Type 1 fiber CSA, urn 2 


8703 ±1858 


81 54 ±2460 


8524±2655 


7756±2268 


Type l/lla fiber CSA, urn 2 


6540 ±2048 


6073 ±2742 


6274±2810 


5838 ±2766 


Type Ila fiber CSA, urn 2 


7107±2580 


7016±2309 


7030±2533 


7001 ±2144 


Type lla/llx fiber CSA, urn 2 


7214±3139 


5760±2795 


6082±3124 


5341 ±2396 


Type llx fiber CSA, um 2 


4694 ±26 13 


3792±2032 


4375±2192 


3257±1799 


CS activity, umol/min/g protein 


63±37 


54±22 


52 ±24 


56±21 


HADH activity, umol/min/g protein 


22±9 


19±7 


19±7 


19±6 


OXPHOS complex 1 subunit content, AU 


2.0 (1.2, 3.9) 


1.5 (0.8, 2.1) 


1.6 (1.0, 1.9) 


1.2 (0.3, 2.4) 


OXPHOS complex II subunit content, AU 


2.8 (1.7, 4.3) 


2.2 (1.5, 3.6) 


2.3 (1.6, 2.5) 


2.2 (1.2, 4.2) 


OXPHOS complex III subunit content, AU 


1 .5 (0.7, 2.0) 


1.2 (0.8, 1.7) 


1 .2 (0.9, 1 .4) 


1 .4 (0.7, 1 .8) 


OXPHOS complex IV subunit content, AU 


4.1 (2.8, 8.1) 


2.4 (1.1, 3.6)* 


2.9 (1.3, 4.1) 


1.9 (0.6, 3.3)' 


OXPHOS complex V subunit content, AU 


0.7 (0.4, 1.1) 


0.5 (0.3, 0.7)* 


0.5 (0.3, 0.8) 


0.6 (0.3, 0.7) 



Values are expressed as mean±SD or median (25th percentile, 75th percentile). Abbreviations: CSA, cross-sectional area; CS, citrate synthase; HADH, 3-hydroxyacyl-CoA 
dehydrogenase; OXPHOS, oxidative phosphorylation; AU, arbitrary units. Significance of difference compared to controls: * p<0.05, ** p<0.01, *** p<0.001. 
doi:1 0.1 371 /journal.pone.00901 50.t002 
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Table 3. Results from the maximal cycle ergometry test. 





Healthy controls 


COPD 


COPDitir. 


COPD D 


Total exercise time, min 


10.0 (8.8, 10.6) 


7.4 (5.5, 9.4)"* 


7.5 (5.5, 8.9)* 


7.0 (5.7, 9.7)* 


Peak work load, Watt 


156 (134, 255) 


72 (56, 96)"* 


76 (66, 100)*** 


66 (53, 96)*** 


Work load, % predicted 


133 (106, 152) 


60 (44, 71)*" 


65 (48, 73)*** 


57 (37, 69)"** 


V'o 2 p eak , ml/min 


1852 (1680, 2785) 


1130 (1025, 1406)*** 


1 162 (1075, 1434)*** 


1084 (901, 1358)""* 


y'°2 P eak. % predicted 


120 (102, 155) 


74 (55, 84)*" 


78 (66, 88)""" 


64 (51, 74)*** 


l/'E.peak, % MVV 


66 (61, 80) 


85 (81, 94)"** 


82 (77, 94)*" 


89 (82, 94)** 


HRp e ak, % HR predicted 


100 (97, 106) 


78 (68, 89)*"* 


77 (66, 89)** 


78 (71, 92)"* 


AWo 2 , ml/min 


1061 (862, 1395) 


736 (607, 873)"* 


776 (601, 875)*" 


727 (605, 889)" 


S p ,02 at W"o 2peak , % 


97 (97, 99) 


92 (87, 95)""" 


94 (93, 96)*" 


87 (85, 90)*"' *** 


AS p ,o 2 , % 


-1.0 (-2.0, 0.0) 


-4.0 (-8.0, -2.3)*"* 


-3.0 (-4.0, 0.0) 


-8.0 (-10.0, -6.5)***' *** 


AS a ,o 2 , % 


0.0 (0.0, 1 .0) 


-1.0 (-4.0, 0.0)* 


-1.0 (-1.8, 1.5) 


-4.0 (-10.0, -3.0)""' ** 



Values are expressed as median (25th percentile, 75th percentile). Abbreviations: V'o 2 , oxygen uptake; I/'e, minute ventilation; MVV, maximal voluntary ventilation, 
calculated as FEV, x37.5; HR, heart rate; HR predicted calculated as (220 — age); AWo 2 , l/'o 2 at anaerobic threshold; S p ,o 2 , oxygen saturation measured via pulse 
oximetry; S a ,o 2 , arterial oxygen saturation. Level of significance of difference vs. controls: *p<0.05, **p<0.01, "*p<0.001. Level of significance of difference COPD D vs. 
COPD ND patients: * # p<0.01, ** # p<0.001. 
doi:1 0.1 371 /journal.pone.00901 50.t003 



The increases in HO-1 and MnSOD gene expression did not 
correlate with exercise-induced oxidative stress-derived products; 
only MDA-protein adducts were slighdy increased (Figure 2C) in 
COPD ND patients. 

Leg muscle desaturation 

Although by definition COPDrj patients had significandy 
reduced end-exercise systemic saturation compared to COPD ND 
patients and controls, no differences were detected in quadriceps 
tissue desaturation or time constant of recovery of muscle oxygen 
saturation among the groups (Table 4). Also, no association was 
found between systemic and local muscle desaturation (Figure 3). 
To assess whether the Oxphen response was influenced by muscle 
oxygen saturation, we tested correlations between Oxphen 
markers and ATSI in the combined group of COPD patients 
and controls, but did not find significant associations. In addition, 
we stratified the COPD patients based on the median value of 
ATSI to high (<-9.1%) and low muscle desaturation (>-9.1%) 
groups, due to lack of a reference value for muscle desaturation. 
The response of Oxphen markers was not different between the 
two groups. 

Discussion 

The mechanisms leading to loss of muscle Oxphen in COPD 
thus far remain unclear. We recently showed that the Oxphen loss 
is not limited to advanced disease but also ongoing in earlier 
phases of the disease, even in absence of muscle wasting. Potential 
determinants such as disease-associated persistent systemic inflam- 
mation and chronic hypoxia are less likely to be causally involved 
in this stage. As exercise is a potent Oxphen stimulant we 
hypothesized that a blunted regulatory Oxphen response to 
exercise would occur in less advanced COPD, resulting from 
either exercise-induced hypoxia or oxidative stress, or both. 
However, we did not demonstrate such a blunted response in 
quadriceps muscle, regardless of co-occurring oxygen desatura- 
tion. Also exercise-induced oxidative stress was not aggravated in 
these patients with less advanced COPD compared to controls. 
Moreover, as the COPD patients exercised at an intensity 
significantly lower than that of healthy controls, it could be 



argued that, adjusted for absolute intensity, the exercise response 
was even larger in the patients as compared to the controls. In line, 
after a constant load test at 80% peak capacity in COPD patients 
and healthy controls, Steiner et al. also did not find significandy 
different exercise-induced changes in absolute magnitude of the 
muscle metabolites ATP, inosine monophosphate, phosphocrea- 
tine and lactate, despite marked differences in workload between 
COPD patients and controls [36]. 

Exercise intensity and response 

One explanation for the comparable exercise responses in the 
COPD patients and healthy controls could be that the induction of 
Oxphen is not linearly related to exercise intensity but rather 
determined by a threshold that has to be exceeded to invoke a 
maximal response. The fact that in our study the anaerobic 
threshold was reached in all healthy controls and almost all COPD 
patients would be supportive of this hypothesis. Congruendy, in 
healthy individuals only exercise above the lactate threshold has 
been shown to significantly induce muscle PGC-lot gene 
expression [37,38]. In COPD patients, however, exercise both 
below and above the lactate threshold increased PGC-lot gene 
expression, although the induction was stronger with the latter 
[38]. Puente-Maestu et al. found increased induction of PGC-la 
mRNA in muscle of COPD patients compared to healthy controls 
after a 45-minutes exercise test at 65% of J /, o 2 peak [38]. Whereas 
the COPD patients in our study had significandy increased PGC- 
loc gene expression after exercise compared to baseline, the 
induction in healthy controls only showed a trend to significance. 
However, we did not observe a significant difference in the 
magnitude of the exercise-associated induction between the 
patients and controls. Puente-Maestu et al. found an association 
of PGC-la mRNA induction with the amount of ROS that was 
estimated to be produced during the exercise test [38]. In our 
study there was no evidence of increased oxidative stress in the 
quadriceps muscle of COPD patients, which together with the lack 
of a differential response of PGC-la mRNA with exercise, is in 
line with the findings by Puente-Maestu. 
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Figure 1. Gene expression responses to exercise are not significantly different between COPD patients and healthy controls, 
regardless of desaturation status. A Gene expression inductions (post/pre) for regulators of Oxphen. PPAR, peroxisome proliferator-activated 
receptor; PGC, PPAR gamma co-activator; PRC, PGC-1 -related co-activator; ERRa, estrogen-related receptor oc; NRF, nuclear respiratory factor; TFAM, 
mitochondrial transcription factor A. B Gene expression inductions for metabolic enzymes involved in citric acid cycle, (3-oxidation and glycolysis. CS, 
citrate synthase; HADH, 3-hydroxyacyl-CoA dehydrogenase; COX4-1, cytochrome c oxidase subunit IV isoform 1; HK-II, hexokinase II; PFK, 
phosphofructokinase muscle isoform. C Gene expression inductions of hypoxia-associated targets. HIF-1oi, hypoxia-inducible factor 1 alpha; VEGF-A, 
vascular endothelial growth factor A; GLUT1, glucose transporter 1. Whiskers indicate minimum and maximum values. Level of significance: p<0.05 
for post vs. pre. 

doi:10.1371/journal.pone.0090150.g001 



Regulatory gene responses to exercise 

Transient induction of PGC- la gene expression by exercise is 
well-described (i.e. [13]). Remarkably, in addition to a mild 
induction of PGC- lot gene expression, we found a strong induction 
of PRC, a PGC-l-related coactivator involved in mitochondrial 
respiratory function. Induction of PRC with exercise has been 
described before in healthy individuals [39], but this is to our 
knowledge the first time that this marker has been measured in 
COPD patients in response to acute exercise. 

Exercise-induced elevation of HIF- 1 a gene expression has been 
proposed to be an important factor in the regulation of adaptive 
gene responses to exercise [40,41]. In addition to HIF-la, we 
found increased expression of some of its target genes, such as 
VEGF-A, HK-II and GLUT1. The inductions of these genes and 
PGC- la and PRC support that there is an ongoing exercise 
response in the muscle. It is thought that induction of the HIF-la 
gene with exercise does not depend on hypoxia per se, but may 
result in increased capillarization and increased glycolytic flux 
within the mitochondria [40]. Inductions of HIF-la gene 
expression have been found to blunt with training [40], and the 
seemingly stronger inductions of HIF- 1 a gene expression and its 
target genes in COPD patients compared to controls could be a 
reflection of the on average lower physical activity level of the 
patients. 



Exercise-associated hypoxia 

In contrast to what we expected, exercise-induced desaturation 
did not blunt the Oxphen response in COPD patients. This could 
indicate that acute local hypoxia does not hamper the Oxphen 
response. Another explanation for the comparable exercise 
responses between the COPD ND and COPD D subgroups is that 
"systemic" oxygen desaturation as measured by finger pulse 
oximetry does not reflect oxygen desaturation locally in the 
exercising limb muscle. To define exercise-induced oxygen 
desaturation, we used the criterion by the American Thoracic 
Society (ATS), which specifies that the drop in saturation during 
exercise should be &4% [26], It has also been suggested that 
desaturation is only clinically relevant when saturation drops 
below 90% [42], which applied to 12 patients in the current study. 
However, the outcomes did not change when we used this 
alternative criterion (results not shown). Indeed, we did not find an 
association between systemic and skeletal muscle desaturation. 
The mismatch between systemic and muscle desaturation could be 
explained by an increased oxygen extraction rate from the 
systemic circulation to maintain oxygen tension in the muscle 
but resulting in decreased arterial oxygen saturation in some 
patients. Nevertheless, we found no significant difference in 
exercise-induced muscle oxygen desaturation levels among the 
groups, including the healthy controls, and no associations of 
muscle desaturation with markers of the Oxphen response. 
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Figure 2. Oxidative stress in quadriceps muscle is only mildly increased in non-desaturating COPD patients; both patients and 
healthy controls show significant induction of antioxidant genes after exercise. A Box plots of the response of heme oxygenase-1 (HO-1) 
and manganese superoxide dismutase (MnSOD) genes to exercise presented as fold induction (post/pre). B Box plot of the exercise-induced fold 
inductions of catalase, Cu/ZnSOD and MnSOD content. C Box plot of the exercise-induced fold inductions of protein carbonylation, tyrosine nitration 
and malondialdehyde (MDA)-protein adducts. Whiskers indicate minimum and maximum values. AU, arbitrary units. Level of significance: p<0.05 for 
post vs. pre. 

doi:1 0.1 371 /journal.pone.00901 50.g002 



Exercise-associated oxidative stress response 

Previous research in advanced COPD demonstrated an 
increased resting and exercise-induced oxidative stress response. 
We also investigated relevant oxidative stress markers in this study 



Table 4. Quadriceps desaturation during cycle ergometry. 
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Values are expressed as median (25th percentile, 75th percentile). 
Abbreviations: TSI, tissue saturation index; x, time constant of recovery of 
muscle oxygen saturation. 
doi:1 0.1 371 /journal.pone.00901 50.t004 
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Figure 3. Change in quadriceps muscle saturation during 
exercise is not associated with systemic desaturation. Reliable 
muscle saturation data was available for 1 1 COPD ND , 9 COPD D patients 
and 1 1 healthy controls. 
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to see if a similar enhancement is seen in less advanced COPD, in 
particular in desaturating patients. No striking differences were 
observed between COPD patients, subgroups and healthy 
controls. In contrast Puente-Maestu et al. recendy described 
increased exercise-induced oxidative stress in mild to moderate 
COPD patients, who, similar to those in our study, had a slight loss 
of muscle Oxphen [43] . However, differences in the intensity of 
exercise (15 minutes constant work rate at ~65% W max vs 
incremental exercise until exhaustion) and the timing of the muscle 
biopsy after the exercise test (1 h vs 4 h) hamper comparison 
between our study and theirs. 

Although we did not find increased oxidative stress-derived 
products after exercise, gene expression levels of the antioxidants 
MnSOD and HO-1 were increased, which is a normal response to 
exercise [28,44]. The induction of MnSOD and HO-1 gene 
expression could indicate that the stimulation of the antioxidant 
system by exercise is not blunted in COPD patients and increased 
formation of reactive oxygen species by exercise is well-balanced 
by the antioxidant system [45]. 

Limitations of the study 

A limitation of this study is that only one biopsy was taken after 
the exercise bout. The timing of the single post biopsy was chosen 
after careful review of the literature focusing on Oxphen gene 
expression after exercise, but no inferences can be made regarding 
differences in the temporal pattern of gene expression changes 
between COPD patients and controls. 

An important constraint of the near-infrared spectroscopy 
(NIRS) technique is the limited optical penetration depth. As such, 
subcutaneous fat disturbs or even obstructs the quadriceps NIRS 
measurement. Therefore, subjects with a large skin fold at the site 
of the NIRS probe were excluded from the NIRS measurement in 
our study. We explicitly report this methodological issue as it often 
seems to be overlooked in the existing literature. 
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